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ABSTRACT 

By combining the 0.12 square degree F814W Hubble Space Telescope (HST) and Spitzer MIPS 24/im 
imaging in the First Look Survey (FLS), we investigate the properties of interacting and merging Mid- 
Infrared bright and faint sources at 0.2 <z< 1.3. We find a marginally significant increase in the pair 
fraction for MIPS 24/im detected, optically selected close pairs, pair fraction^ 0.25 ± 0.10 at z~l, 
in contrast to 0.11 ± 0.08 at z ~0.4, while galaxies below our 24/mi MIPS detection limit show a 
pair fraction consistent with zero at all redshifts. Additionally, 24/im detected galaxies with fluxes 
> 0.1 mjy are on average five times more likely to be in a close galaxy pair between 0.2 <z< 1.3 than 
galaxies below this flux limit. Using the 24/im flux to derive the total Far-IR luminosity we find that 
paired galaxies (early stage mergers) are responsible for 27%±9% of the IR luminosity density resulting 
from star formation at z ~ 1 while morphologically classified (late stage) mergers make up 34% ±11%. 
This implies that 61% ±14% of the infrared luminosity density and in turn ~ 40% of the star formation 
rate density at z ~ 1 can be attributed to galaxies at some stage of a major merger or interaction. We 
argue that, close pairs/mergers in a LIRG/ULIRG phase become increasingly important contributers 
to the IR luminosity and star formation rate density of the Universe at z >0.7. 

Subject headings: galaxies: evolution - galaxies: formation - galaxies: interactions - galaxies: starburst 



1. INTRODUCTION 

Hierarchical models and observations suggest that 
galaxy mergers and interactions play a key role in galaxy 
assembly and star formation, but to what extent is still 
unclear. Studies of gas-r ich mergers in t he local uni- 
verse (e.g., Antennae; see H chweizcr 198 2]) and N- body 
simulations (jMihos Sz He rnauist 19961 lBarnesl 2004) have 
revealed fundamental signatures of the galaxy merger 
process, including tidal tails, multiple nuclei, and vio- 
lent bursts of star formation. While interaction-induced 
star formation is thought to be primarily responsible 
for ultra luminous infrared galaxies (ULIRGs, which 
have Lir > lO 12 !^ ) both locally and a t high redshift 
([Sanders et al.l 119881 : iDasvra et al.l [2006) , luminous in- 
frared galaxies (LIRGs, Lm ~ 10 11 — 10 12 L Q ) appear to 
have multiple driving mechanisms, merger-induced star 
formation being only one. 

Luminous infrared (IR) galaxies are thought to be the 
dominant producers of the cosmic infrared background 
(CIRB), and major contributors to the evolution of the 
cosmic sta r formation rate (CSFR) of galaxies, esp ecially 
at z > 0.7 ([Elbaz et alj|2002t iLe Floc'h et al.ll2005t ). The 
rapid decline from z ~ 1 of the CSFR density has been 
linked to a decline in the merger rate. However, recent 
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close pair studies have suggested that the merger rate has 
remained fairly constant from z ~ 1 ([Bundv et al.l 120041 
iLin et all2004l ). and at z > 0.7 the IR populati on is dom- 
inated by morphologically normal galaxies ([Bell et alJ 
120051: iMelbourne et al.l[200l lLotz et al.ll2006l ). The com- 
bination of these two results suggest that the bulk of star 
formation at z ~ 1 is not driven by major mergers. 

However it must be noted that different merger selec- 
tion criteria probe different stages of the merger pro- 
cess. Quantitati v e meas ure ments of galaxy asymm etry 
([Abraham et all (|1996aH bl): IConselice efail (|2003l )) are 
more likely to probe later stages, while early stage merg- 
ers can be identified by carefully searching for close com- 
panions. There should be some overlap between these 
techniques if galaxy pairs are close enough to have in- 
duced strong tidal interactions, but galaxies in pairs 
could also have normal morphologies, hence if early 
stage mergers are not considered, the impact interac- 
tions/merging have will be underestimated. 

Traditionally, close pa ir studies have been carried out 
in the optical/near-I R (iPatton et all 119971, [20001. 120021: 
Carlberg et "atl 120001: Ite Fevre et al.l 120001 : ILin et all 



20041: iBundv' et al. f i200l ~ However recent investigations 



have begun to explore the Mid-IR properties (star for- 
mation) of galaxy pairs, finding a Mid-IR en hancement 
in pa irs separated by less then ten's of kpe's (|Lin et al.l 
I2006T) . The amount of IR luminosity stemming from in- 
dividual processes (star formation or fueling an AGN) in 
interacting pairs and mergers still remains open. To in- 
vestigate this question we have conducted a study of the 
frequency of MIPS 24/im detected, and undetected close 
optical galaxy pairs and morphologically defined mergers 
in the Spitzer First Look Survey (FLS) 6 . We find that the 

6 For details of the FLS observation plan and the data release, 
see |http://s sc . spit zer . caltech . ed u/fls[ 
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fraction of 24/im detected, optically selected close pairs 
and mergers increases with redshift, and are important 
contributors to the IR luminosity and star formation rate 
density at z ~ 1. 

In the discussion that follows, any calculation re- 
quiring cosmology assumes J7m=0.3, Qa=0.70, and 
H =70kms- 1 Mpc~ 1 . 

2. PHOTOMETRIC AND SPECTROSCOPIC 
OBSERVATIONS 

The Spitzer extragalatic component of the FLS is 
a 3.7 deg 2 region centered around R.A.=17 ,l 18 m 00 s , 
decl.=59°30 00 . Observations of this field were taken 
using all four Infrared Array Camera (IRAC) channels 
(Fazio et al. 2004) and three Multiband Imaging Pho- 
tometer (MIPS) bands (Rieke et al. 2004). Additional 
ground base imag es in u*,g' from CFHT's MegaCam 
(jShim et al.l [20061 ). g', i' data from Palomar 200" LFC 
and NOAO 4-m R and K' band (Fadda et al. 2004; Glass- 
man et al. 2006 in prep) have also been obtained. This 
work focuses on the 0.12 deg 2 ACS-HST F814W imag- 
ing of the verification strip, which has 3er depths in MIPS 
24/im of O.lmJy. Object detection and photometry were 
performed using Sextractor (Bertin & Arnouts 1996). 
Particular care was taken to ensure accurate de-blending 
of galaxies in close proximity to one another, while avoid- 
ing detections of substructure within a single galaxy, con- 
sistent with other reductions of HST imaging with close 
galaxy pairs in mind (Patton et al. 2005). There were 
~59,000 sources extracted within the F814:Wab band 
(hereafter extracted magnitudes referred to as Tab)- Wc 
compared our number counts to those from the Hubble 
Deep Field (HDF) North and South and determined a 
limiting magnitude of Iab ~27.4. 

Using the full MIPS catalog from the FLS we selected 
24/im sources within the area covered by the ACS imag- 
ing (~0.12 deg 2 ). In order to correlate the MIPS ob- 
jects with those identified in the optical we first cross- 
identified sources from the MIPS 24/im sample to the 
IRAC catalog using a tolerance radius of 2.0". This 
choice was primarily motivated by the FWHM of the 
MIPS 24/im (PSF^6") and confirmed by visual inspec- 
tion. Wc then cross-correlated the IRAC/MIPS catalog 
to the ACS sample which we band merged with u*, g' and 
R requiring a positional agreement of <1". When multi- 
ple counterparts were identified, we selected the closest 
object. Ultimately we found 1155 ACS sources also de- 
tected by IRAC and MIPS at 24/im . 

The redshifts used in this study were determined exclu- 
sively from optical spectroscopy. They were obtained by 
cross-correlating the ACS sample, limited to Iab <26.5 
(TV^29,000) with various FLS spectroscopic datasets. 
The vast majority of the included redshifts (>97%) were 
obtained with the Deep Imaging Multi-Object Spectro- 
graph (DEIMOS) on the W.M. Keck II 10-m telescope; 
however, the final sample also included a few redshifts 
based on Sloan Digitized Sky Survey (SDSS) and WIYN 
Hydra/MOS (Marleau et al. 2006 in prep) spectra. 
Galaxies in the FLS Verification region were targeted for 
spectroscopic follow-up during two DEIMOS campaigns 
that bracketed Spitzer's launch. The selection criteria 
for these campaigns are summarized below. For the 
2003 pre-launch campaign, targets were selected based 
on NIR (K s ) and optical (g,R,i) colors. The primary 



sample included sources with K s <20.2, R >19.0 and a 
g, R, i color selection that restricted the numbers of low 
redshift (z<0.6) sources. For the 2004 post-launch cam- 
paign, a purely 24/im selected sample (/24>120 /iJy) was 
targeted for follow-up. The combined Iab distribution 
of targeted and detected sources is shown in Figure Q] 
(top) along with the cumulative redshift identification 
efficiency (bottom). The overall spectroscopic complete- 
ness (defined here as the fraction of targeted sources with 
high quality redshifts) is ~70% for the full sample and 
^80% for sources with Iab <25.0. For a more detailed 
description of the observing strategy, primary selection 
criteria and the o verall flux and redshift distributions see 
iChoi et all (f2006h . 

Since we were exploring the Mid-IR properties of galax- 
ies no optical limit was imposed, instead an IR lumi- 
nosity cut ( L IR > 5.0 x 10 10 or 1.0 x 10 11 L@) was 
used, so that a fair comparison could be made at different 
redshifts. The absolute B magnitude (Mb) distribution 
of the MIPS spectroscopic sample between 0.5<z<1.3 
probes -21<Mb<-19 fairly uniformly and we are not 
strongly biased at higher redshifts. 




'»B ( ma g) 

Fig. 1. — (Upper) The Iab mag distribution for the FLS sample 
(solid), the sample targeted for spectroscopy (long-dash), and those 
where spectroscopic redshifts were acquired (dotted). (Lower) The 
cumulative redshift identification efficiency. The overall spectro- 
scopic completeness is ~ 70% for the full sample and ~80% for 
sources with Iab < 25.0. 

Cross-correlation of the band-merged photometric cat- 
alogs with the redshift samples results in a data set of 
476 sources with Iab < 26.5 and 0.2 <z< 1.3. Of those, 
245 (51%) are MIPS 24/im -detected with a measured 
L IR > 5.0 x 10 10 L Q . The remaining 231 (49%) are non- 
detected at 24/im . Despite the fact that the MIPS and 
non-MIPS galaxies were selected slightly differently, the 
resultant colors of objects with spectroscopic redshifts 
have uniform color properties (Figure [5]). 
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Fig. 2. — Color-color plot where red squares depict the 24/im 
detected objects with spectroscopic redshifts, while blue triangles 
show the undetected 24/im spectroscopic sample. The black dots 
represent the full FLS-ACS catalog for comparison. The 24/im 
detected and undetected spectroscopic samples occupy the same 
color space. Some objects were not detected in all four bands due 
to the field coverage, depths and seeing differences between filters. 



3. STATISTICS OF CLOSE PAIRS, AND MERGERS 

To properly constrain the role interactions and mergers 
play in galaxy evolution, all stages of the process must 
be considered. Typically, merger history analyses utilize 
either pair or structural methods. Galaxies in pairs are 
pre-mergers, or systems undergoing interactions, while 
morphological or structural methods find galaxies that 
have already undergone a merger and are dynamically 
relaxing. 

When discussing the pair fraction, and inferred merger 
rates (as defined in Section 4), it must be noted that these 
measurements are highly dependent on the techniques 
and selection criteria used to identify ongoing mergers, 
especially for galaxies at high redshifts. 

The first method, which we will call the "close pair 
method" , is to count the number of galaxy pairs within 
some projected separation, r pro j CCtc d, and magnitude dif- 
ference (Am). If we assume that these systems will 
merge within a given time-scale due to dynamical fric- 
tion, we can determine the merger rate. Although not 
all pairs will merge, they can potentially trigger star for- 
mation through gravitational interactions. An alterna- 
tive is to select merging systems based o n morphological 
indic ators either by overall appearance (|Le Fevre et alJ 
2000) or computational measurements s uch as asymme- 
try (A), and dumpiness (S) of a system (|Conselice et all 
f2?KfflL 12001 lConselicd[2003li. or Oini coefficient (G), an d 
M 20 parameters ()Abraham et~aT1l2003l : lLotz et al.|[2006h . 
Due to the comparatively limited spatial resolution of 
Spitzer (compared with optical imaging), seaching for 
close galaxy pairs or morphological signatures of interac- 
tion at Mid-IR wavelengths is currently restricted to the 
nearby universe. However, we can correlate optically- 
selected pairs/mergers with global Mid-IR properties and 
investigate the IR activity in these systems out to high 
redshifts. 



3.1. Pair Statistics 

We applied the close pairs technique to identify the av- 
erage number of close companions per galaxy, hereafter 
N c . This measurement is similar in nature to the pair 
fraction when there are infrequent triples or higher order 
N-tuples. Since this is the case here, N c will be occasion- 
ally referred to as the pair fraction. Companions were 
selected using a standard operational close pair defini- 
tion of 5/i _1 kpc< r pro j cc tod < 20/i _1 kpc, and an optical 
magnitude difference (Am) < 1.5 compared to the host 
galaxy, to select nearly equal mass major mergers. The 
term "host" or "primary" galaxy are both used to ref- 
erence the pair member with a measured redshift. The 
inner radius is applied to avoid detection of substructure 
within a galaxy, while the outer 20/i~ 1 kpc limit repre- 
sents the radius within which satellites are expected to 
strongly interact with the halo of the host and merge 
within 0.5-0.9 Gyrs (Patton et al. 1997; Conselice et al 
2003). We find 87 close pairs out of 476 galaxies which 
fulfill these criteria (see Table 1). 

To study the fraction of IR-bright galaxies in pairs, 
we split the pair sample into two sub-sets: those which 
were detected and those undetected with MIPS at 24/im 
down to the flux limits of our survey (0.1 mJy). Figure 
[3] shows a subset of close pairs (both detected and unde- 
tected at 24/im ) with MIPS contours. Due to the small 
separations of close pairs (20/i _1 kpc corresponds to 3.6" 
at z ~ 1) relative to the beam of the MIPS 24/im im- 
ages (FWHM ~6"), there are a few instances (5) where 
only a single 24/im detection is found centered between 
the pair members (see middle left image in Figure [3]) . In 
these cases we assume all 24/im flux is coming from the 
primary galaxy. 

3.2. Field Correction 

Since we have redshift information for only the primary 
galaxy and not the companions we need to consider what 
fraction of these close pairs are a result of random pro- 
jection effects. A field correction was determined using 
two separate methods to account for these close pairs. 
The first assumed the same optical magnitude and red- 
shift distributions independently for both the detected 
and undetected 24/im samples, while the positions were 
randomized. The close pair algorithm was applied to 
50 realizations of these mock catalogs and the average 
N c for each redshift bin was taken to be the pair frac- 
tion expected from random. This assumes the absence 
of clustering. 

We investigated the environments of 24/im detected 
and undetected objects on scales of r pro j e cted ~ 20ft. _1 kpc, 
and found them to be comparable, confirming that the 
increase in pair- fraction of the 24/im detected pairs is not 
because they preferentially lie in clusters. On the other 
hand, there is a weak indication that galaxies detected at 
24/im are more likely to lie in small groups. Since such 
groups may, in some cases, be physical associations, we 
count such cases as separate pairs. However, the number 
of these cases is small, and does not influence our results 
in any significant way. The second method utilizes the 
Iab magnitude distribution of the full photometric cata- 
log (^59,000 sources), and determines the average num- 
ber of companions, within 1.5 mags (Iab), normalized 
to the area covered by 5/i _1 kpc< r pro jcctcd < 20/i -1 kpc. 
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Fig. 3. — A subset of paired galaxies in our sample. Each image 
is 60/i _1 kpc on a side with axes in arcseconds, centered on the 
pair member with the spectroscopic redshift also referred to as the 
primary or host pair member (white circle), while the companion 
is highlighted in a black circle. The two upper rows are close pairs 
which were detected at 24/^m , while the lower set are from the un- 
detected 24/xm paired sample. The 3 — 10cr 24/zm flux contours are 
overlaid. The labels are spectroscopic redshift, and Iab magnitude 
of the primary galaxy. 

The results obtained from the two held correction meth- 
ods agreed within ~ 2%, which is negligible compared to 
the uncertainly in N c . The average of the two methods 
was taken to be the hnal field correction. Both the pair 
catalog and randomly generated catalogs were visually 
inspected for false pairs due to single galaxies being bro- 
ken up into multiple components in the source extraction 
phase, or contaminating stars in the photometric catalog, 
and were removed. 

3.3. Pair Fractions 

The field-corrected optical pair fractions for the 24/im 
detected and undetected sub-samples are presented in 
Figure [4] and Table 1. Errors are computed using 
the jackknife technique (|Efron fc Tibshirani II19861 ). e.g. 
given a sample of N galaxies the variance is given by 
[(N - l)/NJ2i$i] 1/2 - The partial standard deviations, 
Si, are computed for each object by taking the differ- 
ence between N c , the quantity being measured and the 
same quantity with the ith galaxy removed, N C i, such 
that 5i = N c - N ci . 

To allow a more direct comparison to be made between 
the generally lower-luminosity low-z pairs, and those at 
higher redshift, we derived pair fractions for MIPS de- 
tected galaxies with an Lm > 5.0 x 10 10 (approximately 
the IR luminosity of the famous Antenna Galaxies). In 
this way we ensure that the sub-luminous galaxies do not 
strongly influence the pair fractions in the lowest red- 
shift bin. The derived N c for 24/zm detected close pairs 
is - 11% ± 8% atz~ 0.4 and increases to 25% ± 10% 
at z~l. In contrast, close pairs with no 24 ^m detection 
show no increase with redshift and have pair fractions 
consistent with zero at all redshifts. The higher pair 
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Fig. 4. — The field corrected pair fraction N c as a function of 
redshift as measured in the optical and IR. The stars represent the 
measurement from our 24/^m detected sample, triangles depict the 
field corrected pair fraction of the undetected 24/xm sub-set, and 
filled circles show the combined pair fraction of the two samples. 
Other optically-determined pair fractions appear as open squares 
(DEEP2 Fie ld 1) and triangles (Field 2) IILin et al J 120041, cross 
for CNOC2 IPatton et aUI2002D . while the dashed line shows their 
best power law fit of (l + z) m (to = 1 08 ± 0.40). The near-IR pair 
fraction determined by Bundy et al. (2004) is shown by diamonds. 
Errors for this work are derived using jacknife statistics, while the 
IR pair fraction errors implemented counting statistics and DEEP2, 
CNOC2 errors are determined via bootstrap. A Lju > 5.0 X 10 10 
limit was imposed on the 24/xm detected close pairs. Note that 
each work imposes a slightly different luminosity and stellar mass 
limit. 

fraction of MIPS bright sources is marginally significant 
due to the small number of sources in the highest red- 
shift bin, more MIR selected samples between z =1-1.5 
are required to strengthen our findings. 

We would like to be able to rule out the possibility 
that N c is biased by the brightest IR sources at z > 1, 
since merger fra ctions change as a function of lum i- 
nosity and mass (|Conselice et al.ll2003 Ku et al.ll200l . 
To address this we placed a higher IR luminosity limit 
{Lib, > 7.0 x 10 11 ) on the sample, so that at z > 0.7 
the same populations were being probed (optically we 
are probing -22 < Mb ~ -19). We still find an increase in 
N c from the lower (0.8 < z < 1.0) to the higher (z > 1) 
redshift bins of similar magnitude compared to when the 
lower IR limit (Lj_r > 5.0 x 10 10 ) was used. Therefore 
the increase in N c found at z > 1 is likely not a result of 
merely probing brighter IR systems but rather due to a 
physical increase in the merge rate for the 24/Ltm popula- 
tion, however deeper 24/im imaging and spectroscopy are 
required to confirm this. When we consider the averaged 
pair fraction over 0.2 < z < 1.3 for the 24/xm detected 
sample we find that galaxies above a flux limit of 0.1 mJy 
are five times more likely to be in a close galaxy pair, than 
those below this limit. 

3.4. Morphological Mergers 

To explore the structural components of galaxies in our 
sample we used the CAS (Concentration, A symmetry, 
dumpiness) quantitative cl a ssification system (iConselicd 
fl997t iConselice et alJl2000l [200l iConselicd [200311 . and 
visual classifications. To measure the merger fraction us- 



Interacting & Merging Galaxies: Spitzer's View 5 



TABLE 1 
FLS Close Pair Statistics 



z 


Ngal 






N c 


K 


24/im Detected 


0.2-0.5 


32 


0.188 (6) 


0.078 (2.5) 


0.110 ±0.083 


0.83 


0.5-0.80 


82 


0.171 (14) 


0.057 (4.7) 


0.114 ±0.040 


0.93 


0.80-1.0 


82 


0.122 (10) 


0.029 (2.4) 


0.093 ±0.038 


0.90 


1.0-1.3 


49 


0.429 (21) 


0.182 (8.9) 


0.247 ±0.086 


0.67 


24/im Undetected 


0.2-0.5 


44 


0.136 (6) 


0.102 (4.5) 


0.034±0.052 


1.00 


0.5-0.80 


76 


0.132 (10) 


0.134 (10.2) 


0±0.039 


1.00 


0.80-1.0 


56 


0.214 (12) 


0.193 (10.8) 


0.021±0.064 


0.83 


1.0-1.3 


55 


0.145 (8) 


0.180 (9.9) 


±0.065 


0.75 



Note. — N ga i is the number of galaxies with a spectroscopic redshift, 
is the number of companions per host fulfilling our pair criteria, while is 
the number of projected companions per host from the field. The corrected 
fraction of companions per host is given as, (N -N^}) / N ga i , with errors being 
determined using a jacknifc technique. Numbers appearing in parentheses refers 
to the number of close pairs in the respective redshift bins. Undetected at 24 )im 
refers to sources below the limits of our survey (0.1 mjy). The constant k is 
the fractional number of mergers per host galaxy. A Ljr > 5.0 X 10 10 limit 
was imposed. 



ing structural classifications we visually inspected the full 
24/im detected spectroscopic catalog with the following 
groupings: early type (E, SO), mid- types (Sa-Sb), late- 
types (Sc-irr), compact systems, disturbed disks, and 
mergers. The methodology for carrying out th i s clas - 
sification is described in detail in lConselice et all (2005). 
Basically, each galaxy was viewed on a computer screen 
and classified into one of our types. Overall we find that 
55% ± 5% of 24/ mi detected galax i es are disks , whic h 
is consistent with iBell et al.l (|2005l ): lLotz et all (|2006l ). 
while 26% ±5% are merging systems and ~ 6% were clas- 
sified as disturbed disks and are possible minor mergers. 
A fraction of the disk-dominated objects do show some 
visual signs of a morphological disturbance, or are in a 
pair, as we will discuss later in this paper. 

Galaxies undergoing a major merger event can also 
generally be identi fied by their large asymmetr ies in the 
rest frame optical ([Conselice et all 120001 [20031) . We de- 
fined a major merger as a galaxy having an asymmetry 
(A) > 0.35 and Iab < 26.5 (see Figure [5] for exam- 
ples). This limit has been shown to be a clean way to 
find galaxy mergers, withou t significant con tamination 
from non-merging galaxies (|Conselicd [2003]) . Figure M 
shows how the merger fraction for CAS defined mergers 
evolves as a function of redshift for both 24/im detected 
objects (top panel) and LIRG/ULIRG galaxies (bottom 
panel). As with the 24/im detected close pair sample 
there i s an elevated merger fraction compa red to other 
works (jCassata et al.|[2005l lLotz et aT1l2006l ) in which no 
24/im flux limit was imposed, and a slight indication of 
evolution with redshift, but it is statistically consistent 
with 77i ~ 1.0 (dashed line), where m is the slope of a 
power-law of form (1 + z) rn later discussed in §5. 

We also performed a CAS analysis of our close pairs 
sample which revealed that 24/im detected pairs are no- 
tably more asymmetric than the undetected-MIPS close 
pairs (Figure [7]), suggesting that interactions and colli- 
sions may play a role in their IR activity. If the 24/im 
detected close pairs were generally of a different mor- 
phological classification than those pairs undetected at 
24/im the discrepancy in the asymmetries could be ex- 
plained. To address this issue each close pair was visually 




in arcseconds. The 3 — 10o- 24 flux contours are overlaid. The 
labels include spectroscopic redshift, and 7^81411^43 magnitude. 

inspected and classified by four of the authors to be ei- 
ther disk or bulge-dominated. We find that 81% of the 
24/xm pairs have disk morphologies while 74% of the un- 
detected 24/im hosts were also disk dominated, hence the 
discrepancy between the asymmetries of the two groups 
is not caused by classification differences, but rather is a 
physical effect. 

4. MERGER RATES 

One of the goals of studying mergers and interactions 
is to determine how the galaxy merger rate evolves with 
redshift. Most studies of galaxy mergers involve deter- 
mining the merger fraction, yet the merger rate, which is 
defined as the number of galaxies merging per unit time 
per unit volume, is a more physical quantity that can 
be used to determine the full merger history. The rate 
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Fig. 6. — Merger fraction as a function of redshift, using quan- 
titative morphological criteria. The filled squares represent mea- 
surements from FLS (this work) of sources with a_24/im detec- 
tion. The "x's" mark the resu lts of Conselice et al. (200 31), triangles 
depict [Cassata et al. (2005), and diamonds represent ILotz et al.l 
(2006). Top panel shows the merger fraction of other studies with 
no 24/^m criteria imposed, while the measurements from this work 
are mergers with Lir > 5.0 X 10 10 Lq. The dashed lines shows the 
best fit of (l + z) m using all points (from the other studies) with no 
MIPS limit imposed (m = 1.08). Bottom panel shows the merger 
fraction for LIRG/ULIRG galaxies (L IR > 1.0 X 10 n L Q ). Error 
bars derived using Poisson statistics. 
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Fig. 7. — Asymmetry - Concentration diagram for the 24/mi de- 
tected close pairs (squares) and undetected pairs (triangles). The 
long dashed (vertical) line separates merging and non-merging sys- 
tems (A> 0.35 is considered a merger), while the dash-dotted lines 
se parate early (upp er) to mid to late (lower) type galaxies defined 
bv lCon sclicc (2003). Generally close pairs are mid/late type galax- 
ies, and some would also be morphologically classified as a merger. 

in which galaxies merge also affects the mass function of 
galaxies, and is likely linked to the cosmic star forma- 
tion rate. Since we are considering a very broad range in 
the merger process, from early-stage or pre-mergers se- 
lected via close galaxy pairs, and later-stage mergers cho- 
sen based on morphological criteria, we must be careful 
when determining their respective merger rates, as the 
time-scales for these processes are all different. 

There are two variations of the merger rate definition. 
The first is the number of mergers that a galaxy will 
undergo per unit time (K ms ), and the second is the total 



number of mergers taking place per unit time per unit co- 
moving volume (Wlmgv)- Since we are primarily interested 
in mergers which are also Mid-IR bright systems we will 
have to restrict ourselves to measuring 5R mg because the 
evolution of the 24/im luminosity function with redshift 
is currently not well constrained, and our redshifts are 
not complete enough to reconstruct this evolution. 

In order to determine JR mg we need to identify sys- 
tems which are destined to merge. We have approached 
this measurement from three different perspectives, close 
pairs to select pre-mergers or interactions, visual inspec- 
tion to select interactions after the first passage, and 
late stage mergers, as well as CAS criteria which quan- 
titatively selects for later stage mergers. By combin- 
ing information about the number of ongoing mergers 
(N m ) and the time-scales, (T mg ) on which they will un- 
dergo said merger, one can estimate an overall merger 
rate 3i OTff = N m /T mg . Each method of identifying merg- 
ers/interactions is capturing a different snapshot of the 
merger process, each with different merger timescales. 

The value of N c is directly proportional to the num- 
ber of mergers per galaxy (N m ), such that N m = nN c 
(re is a constant relating to the number of mergers per 
galaxy). Hence, the merger rate detemined using close 
galaxy pairs is given by 5R mfl = nN c /T mg . The value 
of re depends on the nature of the merging systems un- 
der consideration. If one were to identify a pure set of 
galaxy pairs each consisting of one companion undergo- 
ing a merger, then n — 1.0. In our case it exclusively 
accounts for close pairs which are in doubles and per- 
haps higher order N- tuples. Our definiti on of k differs 
by a factor two from iPatton et alj (|2000f) which in this 
instance would have re = 0.5 since they have redshifts 
for both pair members and one merger is made up of 
two companions. We have redshift information for only 
one pair member, therefore one merger is made up of a 
primary and one companion. The merger rate equation 
for merging galaxies selected by visual classification and 
CAS parameters is simply 5R„ lg = f mg /T mg , where f gm 
is the galaxy merger fraction. 

Before the merger fraction can be used to calculate 
the merger rate we need to understand the time-scale 
in which a merger occurs. Each technique of identifying 
mergers has a different time-scale since each is sensitive 
to a different interval of the merger process. There are 
two main methods that have been used to estimate the 
time-scale of a merger: dynamical friction arguments, 
and N-body models. The details of t hese methods are be - 
yon d the scope of thi s paper but see lPatton et alj (2000) 
and iConselicd (|2006l ) for a review. We take the aver- 
age merger time-scale for a set of close companions of 
roughly equal mass to mer ge as ~0.5 Gyr±0.25, derived 
from d ynamical argume nts (jPatton et alj|2000t IConselicel 
120061 ). IConselicd (|2006D showed through N-body simula- 
tions that visual classification selects on-going mergers 
over a longer time-scale (1.0 Gyr±0.25) since the hu- 
man eye detects both early and later stage mergers, while 
the asymmetry of a galaxy is sensitive to 0.41 Gyr±0.17 
(jConselicdl2006l) of the merger sequence. 
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5. THE EVOLUTION OF THE GALAXY MERGER RATE 
0-2 < Z < 1.3 

Within the past two decades numerous studies have 
been performed to estimate the evolution of the galaxy 
merger fract i on, u s ing both the close pair technique 
(iZepf fc Kool [1 989: Bu rkev et all 119 94; Carlbcr g~etaT1 
1994 lYee fc Ellingsonl 119951; IPatton et alj [l997[ 120001 ; 
Le Fevre et al.ll200Ct iLin et al.ll2004t IBundv et al.ll2004h 
and morpholog i cal p a rameters (iLe Fevre et al.l 120001; 
Conselice et~aTT 120031 : lLaverv et all 12004 iLotzetall 
20061) . Evolution in the galaxy merger fraction is of- 
ten parameterized by a power-law of form (1 + z) m , and 
has yielded a wide range of results, spanning 0<to<5. 
The large spread in values is in part due to the different 
selection criteria used to identify merging systems and 
biases from optical c ontam ination or redshift complete- 
ness. IPatton et al.l (|1997| ) considered these biases and 
demonstrated that most results to that date were con- 
sistent with their estimate of m = 2.9 ±0.9. Recently , 
optical and near-I R close pair studies (|Lin et alJ l2004t 
IBundv et al.l l2004h have derived merger fractions with 
little redshift evolution (m~l), a s have some morp holog- 
ical studies using (G), and M 2 o (|Lotz et al.ll2006|) . 

When we consider all the close pairs identified in our 
sample, both those detected at 24/xm and not, we find a 
merger fraction and rate consistent with recent studies 
showing little redshift evolution. However, when we sep- 
arate the pair sample into systems with a 24/im detection 
above 0.1 mJy, and those below it, we do see a stronger 
evolution of N c with redshift, (recall that N c oc 5R mff ) and 
therefore also in the merger fraction and rate (Figure [8|). 
Similarly, visually classified mergers and those identified 
via asymmetry levels (A > 0.35) using the CAS param- 
eters, also show redshift evolution in the merger fraction 
and rate. The merger fraction computed using the differ- 
ent methods are in good agreement when normalized by 
their respective time-scales, reinforcing the idea that we 
are probing different phases of the merger process. Con- 
sidering all three merger selection techniques we find the 
best fit of the merger rate parameterized by K(0)(l + z) m 
to be 0.077±0.045, 2.12 ±0.93, with a reduced x 2 =0.39. 
This result suggests that when one considers a sample of 
close galaxy pairs solely on their optical fluxes, brighter 
than Mb =-19, little evolution of the merger rate with 
redshift is found. However close pairs emitting 24/im 
flux exhibit an increase in the merger rate with redshift. 
The infrared luminosity limit (Lib > 10 11 !/©) imposed 
on the close pairs and mergers allows us to primarily 
probe systems in a LIRG/ULIRG phase at z > 0.4 (see 
next section for details). The increase of the merger frac- 
tion and rate of this population of galaxies coupled with 
the fact that LIRG /ULIRG galaxies dom inate the SFR 
density at z > 0.7 (|Le Floc'h et alj|2005h suggests that 
merging does in fact play an increasingly important role 
in star formation out to z ~ 1. 

6. TOTAL INFRARED LUMINOSITIES OF MERGERS 

One way to quantify the role merging galaxies play in 
triggering star formation is to investigate their contribu- 
tion to IR luminosity densities. Infrared luminosities (8- 
1000 fim ) were calculated utilizi ng the 24/im fluxes an d 
two different tem plate methods: ICharv fc Elbazl 1 2001): 
iDale et al.l (l200l in a similar manner as iLe Floc' h et al.l 
( 20051) for bhe full MIPS 24/im spectroscopic sample (Fig- 
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Fig. 8. — The number of mergers per galaxy (Lm > 5x 10 10 ) as a 
function of redshift. Three merger /interaction selection techniques 
are applied, close pairs (squares), CAS criteria (stars), and visual 
classification (triangles), while merger rates for the combined (w/ 
and w/o 24^tm detections) using the close pairs method are shown 
with circles. The long dashed curve is the best fit of the form 
(1 + z) m , using the FLS data for the three techniques; the dot- 
dashed curve represents the best fit for the combined total close 
pairs (MIPS and non-MIPS pairs). 

ure[5]). MIPS pairs and mergers share a similar luminos- 
ity distribution to 24/im bright field galaxies, although 
red-AGN seem generally mor e luminous which is in part 
due to template mismatches (|Chary fc Elhadl200l . 

The Lm of a galaxy is a combined measure of the 
reprocessed UV photons intercepted by dust from mas- 
sive young stars and AGN. Therefore to investigate the 
contribution an interacting or merging galaxy makes to- 
wards the total Ljf> density from star formation alone 
we must first remove AGN from our sample. Due to 
the nonuniform rest-frame spectral coverage of our sam- 
ple we rely on the f our-band IRAC color selection used 
by (jLacv et al.ll2004[ ) to identify and remove AGN can- 
didates (Figure [TOj) . Over the modest redshift range of 
our sample, this method is still effective at separating 
IR-warm AGN from starburst systems. We find an AGN 
contamination rate of ~ 12% for the full 24/im sample, 
while ~ 14% of the hosts in a pair or merger were char- 
acterized as AGN. 

With AGN candidate objects removed we can infer the 
contribution to the Lib. density from star formation com- 
ing from 24/im galaxies in a interaction/merger as a func- 
tion of redshift. We derive the number of statistically 
"real" galaxy pairs from our pair fraction result at each 
redshift interval and determine the total Ljb density 
from close pairs which is in turn divided by the Lib. den- 
sity from the whole sample. We find that paired galax- 
ies (Lib > 10 n £©) are responsible for 27%+^ of the 
IR background stemming from star formation at z ~ 1. 
Since we only know the redshift of the host galaxy we 
select "real" close pairs in a statistical sense, and derive 
error bars for the close pairs contribution by the spread of 
50 realizations of the Lib. density from different combina- 
tions of 24/im galaxy pairs. We also applied this analysis 
to CAS and visually classified mergers, which make up an 
additional ~12%, and ^22% of the IR luminosity density 
respectively. Naturally there is a some overlap in mergers 
identified through close pair criteria and morphological 
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Fig. 9. — Infrared luminosity L IR yg > _ iooo/jm] vs - redshift for 
the MIPS 24 (im sample with spectroscopic redshifts, broken down 
into galaxies in a close pair (red triangle), CAS mergers (open 
green diamonds), AGN candidates (yellow star), AGN candidates 
in a close pair (blue square), AGN mergers (open black square), 
and field galaxies (black circles). LIRGs lie above the horizontal 
dashed line at Ljjj > 10 11 Lq 

parameters, since interacting pairs can exhibit tidal tails 
and asymmetric structures, causing them to also be iden- 
tified morphologically as mergers. We found that 37% of 
CAS defined mergers were also in a close pair, and 31% 
of visually identified mergers were also classified by CAS 
as merging. In cases where a merging system was identi- 
fied using multiple techniques it's contribution was only 
counted once. For example, if a merger identified mor- 
phologically (either through CAS or visual inspection) is 
also in a close pair it is removed from the morphologi- 
cal merger catalog, or if a CAS merger is also identified 
visually the merger is removed from the visual merger 
catalog. This insures that no close pair or merger is 
counted more than once when deriving the contribution 
from interactions and mergers to the IR luminosity den- 
sity. 
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Fig. 10. — IRAC color-color plot using the region in the FLS cov- 
ered by the ACS imaging. Circles represent objects with spectro- 
scopic redshifts and detections in all four IRAC channels. The red 
triangles indicate objects which have met t he color criteria ( shown 
by the dashed line) of an AGN candidate IILacv et al. 2004)- The 
green stars and blue squares depict objects in a close galaxy pair 
or merger whose host was also flagged as an AGN candidate. 
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Fig. 11. — The fraction of Lm density (as a result of star forma- 
tion) as a function of redshift coming from LIRG/ULIRG galaxies 
in a close pair (pre-merger phase) shown by triangles, or more ad- 
vanced stage mergers defined morphologically (circles). The star 
symbol indicates the total combined contribution from close pairs 
and CAS mergers, while squares depict the total from close pairs, 
CAS, and visually classified mergers. Note an infrared limit of 
Lir > 1O 11 L0 was imposed. 

The combination of these three merger selection tech- 
niques identifies a large range in the merger process, from 
pre-merger to late stage mergers, implying that ^60% 
of the infrared luminosity density at 2~1 can be at- 
tributed to galaxies involved in some stage of a major 
merger fFigure fTTj) . The remaining ^40% of the IR back- 
ground from LIRG/ULIRGs is likely to predominately 
come from active, isolated gas-rich star-forming spirals, 
with some contribution from minor mergers. 

If we exclude visually classified mergers the close 
pair/merger contribut ion to the IR de nsity is ~38%, in 
good agreement with iLin et all (|2006l ) who estimate a 
moderate contribution from interacting and merging sys- 
tems of < 36%. It must be noted however that neither 
ILin et all (|2006[ ) or this work have considered the contri- 
bution from minor mergers and are therefore lower limits. 



6.1. Star Formation in Mergers & Interactions 

An important and highly debated question is: how im- 
portant are galaxy mergers in understanding the dra- 
matic decline of the cosmic SFR density from z ~ 
1 to the present day? It has been well established 
that mergers and i nteractions can induce violent bursts 
of star formation (Schwcizer 1982; Bar ton et al.l l2000l : 
iMihos fc Hernquisdll99d lCox et al.ll2006h . So to investi- 
gate this contribution we derived the SFR for our 24^m 
detected close pairs and mergers, using their Lj^. The 
infrared luminosity of a galaxy is a star formation rate 
tracer which is unaffected by the extinction of dust. The 
dominant heat sources of most dusty, high-opacity sys- 
tems such as LIRGs and starbursts is stellar radiation 
from young stars. In these types of systems the Lj_r 
can be co n verted into a SFR using the calibration of 
iKennicuttl lfl998h . SFR IR = 4.5 x 10^ 44 L /fl (ergs s" 1 ), 
where Ljr is the integrated luminosity from 8-1000/im 
as determined in section 6.0. 

We estimated the contribution mergers and interac- 
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tions above Lm > 10 11 L© make to the SFR density at 
z ~1 in two ways. The first is simply to consider their 
contribution to the Lir density which is a star formation 
tracer. Section 6.0 determined that mergers and interac- 
tions at z~l (above Lib, > lO 11 ^©) are responsible for 
40-60% of the IR lumi nosity density. Using the results of 
iLe Floc'h et~aTI (j2005ft which showed that z >0.7 LIRGs 
produce ~ 70% of the star formation rate density, we 
can infer that mergers and interactions in LIRG/ULIRG 
phase would be responsible for ~ 30 — 40% (0.6 x 70%) of 
the SFR density at z ~1, since IR activity traces dusty 
star formation. 

The second more detailed approach utilizes the SFR 
density directly arising from our sample of mergers and 
interactions. At 1.0<2<1.3 we find that 59% (12 close 
pairs, and 17 later stage mergers) of galaxies detected 
at 24/zm are involved in some stage of an interaction or 
merger. In this same redshift range our sample is insensi- 
tive to galaxies with IR luminosities <10 115 . To correct 
for this we derived a scaling factor (~7) simply by com- 
paring the number of observed objects of a given Lm 
in a spe cific redshift range t o the number expected from 
models (Laga che et al.ll2004[ ) . However to go any further 
we must assume that our spectroscopic sample is repre- 
sentative of this population at ~1, and by all accounts 
this appears to be true. Using the derived pair fraction 
we can then infer the total number of major mergers and 
interactions occurring (fulfilling our criteria) in a given 
volume and Lir limit. The lower limit of the SFR den- 
sity at z ~1 from merging and interacting galaxies is 
found to be 0.066 M0yr _1 Mpc _1 . Using the extinction 
correct ed "Lilly-Madau" plot (0.1585 Moyr^Mpc^ 1 at 
z =1) (jThompson et al.ll2001ft we find that mergers and 
interactions are responsible for at least 42% of the SFR 
density at z ~1 (assuming mergers contribute 60% of the 
IR density). Both approaches are in good agreement, and 
are only a lower limit, since objects flagged as AGN were 
not considered even though some of their Lm is a result 
of star formation, and minor mergers which have been 
shown to also induce bursts of star formation were not 
included. 

These results have interesting implications for the 
physical mechanisms that drive the decline in the cos- 
mic SFR (CSFR) density from z ~ 1 to present day. 
They suggest that when all stages of the merger process 
are considered (pre-merger to later stage merger) ma- 
jor interactions and mergers contribute close to half of 
the z ~ 1 SFR density, and the decline in the num- 
ber of 24/xm detected mergers/inteactions is a signifi- 
cant, but perhaps not the primary driver for the decline 
in the cosm ic SFR. This con c lusion differs in i n terpre - 
tation from iBell et al.l (120051 ) : iMelbqurne et aTl 12005); 
IWolf et alJ (12005ft ; iLin et all (|2006ft : lLotz et alj (|2006n , 
which generally suggest that the evolution of the merger 
rate is not a significant underlying cause of the decline 
in the cosmic SFR, but rather a strong decrease in the 
SFR of morphologically undisturbed spiral galaxies is the 
dominant mechanism. Their results do not preclude the 
possibility that their "star forming (undisturbed) disks" 
could be in widely separated pairs, and when we only 
consider quantitatively defined morphological mergers 
our results are consistent with theirs stressing the impor- 
tance of considering the merger process in its entirety. It 



must also be mentioned that we are probing to higher 
redshifts than lBell et~afl (|2005l ). which found that major 
galaxy mergers account for <30% of the IR luminosity 
density at z~0.7, consistent with our findings of 35% 
at that redshift. Our results also agree that at z ~0.7 
isolated undisturbed spiral galaxies are a primary con- 
tributor, however, the influence shifts to interactions and 
mergers at z>0.7. 

Our findings point to an increased importance of MIPS 
bright interactions and mergers to the IR luminosity den- 
sity and SFR density at z>0.7. This conclusion is not 
hampered by the small statistics of the z >1 bin. Fig- 
ure l 111 shows the IR luminosity density contribution from 
interactions/mergers at z~0.7 to be ^37% and 52% at 
z^0.9, reinforcing this increasing trend. 

7. DISCUSSION 

Using a spectroscopic sample of field galaxies from the 
ACS component of the FLS and dividing it into two sub- 
sets, those with a 24/im detection (above 0.1 mJy) and 
those without (or below) we identified optically merg- 
ing/interacting systems via close pair statistics and mor- 
phological methods. We find that roughly 25% of galax- 
ies emitting at 24/xm have a close companion at z ~1 
while at z ~0.5 only ~ 11% are in pairs. In contrast, 
those undetected at MIPS 24/xm showed a pair fraction 
consistent with zero at all redshifts (0.2<z<1.3). On 
average MIPS 24/im galaxies are five times more likely 
than non-MIPS sources to have a close companion over 
0.2<z<1.3. When the samples are combined (regardless 
of 24/im flux ) we find pair fractions consistent w ith previ- 
ous studies (|Lin et al.ll200l iBundv et aI1l2004ft showing 
little evolution with redshift. 

An important and open question is the cause of star 
formation in LIRG galaxies at high-z. Some morpho- 
logical studies have suggested that since at least half 
of th e LIRG ga l axies exhibit disk domi nated morpholo- 
gies (|Bell et alJ 120051: lLotz et al.1 120061) at z >0.7 and 
low non-evolving merger fractions (|Lotz et al.ll2006l ) that 
the driver of IR activity in high-z LIRGs is from on- 
going star-formation from isolated gas-rich spirals and 
not merger or interaction induced. One bias of HST mor- 
phological studies involving the identification of merg- 
ing/intereacting systems is the limitation of detecting low 
surface brightness features such as tidal tails caused by 
close interactions, which can lead to an underestimate of 
the importance of mergers in the evolution of galaxies 
at z <1. Ultimately both close pair and morphological 
techniques must be applied and considered, to obtain a 
complete major merger timeline. Our analysis is the first 
to probe merger rate evolution combining close pairs and 
later stage mergers while considering the IR activity of 
these systems. 

We find that close pair statistics, visually classified 
mergers, and those identified via quantitative CAS pa- 
rameters all showed similar evolution in their merger 
rates. Fitting the merger rate evolution function Sft(z) oc 
(1 + z) rn for 24/im detected mergers above 0.1 mJy, we 
find m = 2.12 ± 0.93. This result agrees with previous 
claims of an increase (m > 2) of the merger rate out 
to z ~ 1 dPatton et al.lll997l I2000L ILe Fevre et alJ l200d: 
IConselice et alj|2003t ICassata et al.ll2005D . However this 
evolution is not seen when IR faint (< 0.1 mJy) merg- 
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ers are included, suggesting that it is the LIRG-merger 
population that is evolving with redshift. 

The Mid-IR emission of LIRGs is indicative of dust en- 
shrouded star formation (and some AGN activity), and 
at z > 0.7 they dominate the IR luminosity density and 
in turn the volume-averaged star formation rate density 
at z ~ 1. We estimate that close galaxy pairs are respon- 
sible for ~ 27% of the IR luminosity density resulting 
from star formation at z ~ 1, while later stage mergers 
contribute ~ 35%. This implies that 40-60% of the in- 
frared luminosity density at z ~ 1 can be attributed to 
galaxies involved in some stage of a major merger, indi- 
cating that merger-driven star formation is responsible 
for 30-40% of the star formation density at z ~ 1. This 
value is a lower limit since minor mergers and interac- 
tions/mergers with an AGN were not considered. 

Ultimately, our findings suggest that interactions and 
mergers of LIRG phase galaxies play an increasingly im- 
portant role in both the IR luminosity and SFR density 
from z > 0.7 out to z ~ 1.3, and are vital to our under- 



standing of the evolution and mass assembly of luminous 
IR galaxies. 
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